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Abstract 

We consider here probabilistic models of transportation flows. The main goal of this introduction is rather 
not to present various techniques for problem solving but to present some intuition to invent adequate and 
natural models having visual simplicity and simple (but rigorous) formulation, the main objects being cars 
not abstract flows. The papers consists of three parts. First part considers mainly linear flows on short 
time scale - dynamics of the flow due to changing driver behavior. Second part studies linear flow on longer 
time scales - individual car trajectory from entry to exit from the road. Part three considers collective car 
movement in complex transport networks. 
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Mathematical models of car traffic can be very different: from partial differential equations to modern com- 
puter graphics where points move on video along the edges of some graph. We consider here only probabilistic 
models. The main goal of this introduction is rather not to present various techniques for problem solving but 
to present some intuition to invent adequate and natural models having visual simplicity and simple (but rig- 
orous) formulation, the main objects being cars not abstract flows. Moreover, the postulates and parameters in 
these models should allow statistical verification and estimation, and should not use doubtful physical analogies. 
Probabilistic models should be tightly related to psychology of drivers, if the drivers are not robots. There is 
no such theory, and we present an attempt to start it. 

Probabilistic approach to transportation problems exists already more than 50 years, see [TJISE], however 
in this paper we follow newer probability theories and consider more difficult problems. There are however some 
works which, for various reasons, we did not include in our introduction, see for example [121 114) 

The papers consists of three parts. First part considers mainly linear flows on short time scale - dynamics of 
the flow due to changing driver behavior. Second part studies linear flow on longer time scales - individual car 
trajectory from entry to exit from the road. Part three considers collective car movement in complex transport 
networks. 



1 Car flows 



1.1 Marked point fields 

Current is normally the mean number J of cars crossing, per unit of time and in the given direction, some 
section of the road. Flow is normally either a static random configuration 

... < Xi < Xi-\ < ... 

of cars at a given time moment, or a probability measure on the set of car trajectories {xi(t)}. 

What is a car configuration In this paper the maximally detailed description for the configuration of cars 
at a given time moment will be as follows. The car is specified by some index a. For example, a road (route) 
with k (traffic) lanes 1,2, ...,fc , is represented by k lines, parallel to x-axis. Then the index a — (m,i) marks 
out the i-th car on the lane m. Index i enumerates cars on the lane so that i-th car follows car with index i — 
Let d a the length of the car indexed by a, x a (t) - its coordinate (for example, of the front buffer). The cars 
move in the positive x-direction. Further on we omit the lane index (the reader can add it whenever necessary) 
and use only index i. 
Denote 

df(t) = x l -i{t) - Xi(t) - di-i 
the distance of the car i to the previous car at time t. The distance of the car to the following car 

dr(t)=df +1 (t) 

is also important for the driver. 

How probabilities on the car configurations are introduced Formally, the point flow on the real line 
R is given by the probability measure on the set of all countable locally finite (that is finite on any bounded 
interval) subsets of R. Otherwise speaking, it is given by compatible systems of probabilities 

P(I\ , k\ ; I n , k n ) 

of the events that in the intervals Ij,j — 1, ...,n there are exactly kj particles. 

The main question of course is how to specify this system more concretely. There are two big parts of 
probability theory which use different ways of doing this. These are renewal theory (see for example, [?>]) and 
Gibbs point fields theory [23, 24J. The first is essentially simpler but works only in one-dimensional case. The 
second has deep relations with physics, is suitable also for multidimensional point configurations but is more 
complicated, and we shall not touch it here. 

The simplest point flow is the Poisson flow, see [30]. Most natural way to understand it is as follows. 
Consider the interval [-N, N] and throw (choose) on it independently and randomly (more exactly, uniformly) 
M = [pN] points, where p > is some constant called density. It is easy to calculate the binomial probability 
Pn,m(^j I) that exactly k points will be inside the finite interval /. This probability tends, as N — > oo, to the 
Poisson expression (Poisson distribution for given I) 

P(fc,J) = M^e-'l'l. 
More general flows can be easily constructed on the half line [0, oo). Namely, random points 

xq = 0, X\) x n , ... 

are defined as the sums 

of i.i.d. (independent identically distributed) random variables £j > 0, i = 1, 2, with the distribution function 
G(x) = Pr(ti < x). 

To define a translation invariant point flow on all real line one problem remains - how to choose point 0, 
from which one can consecutively place independent variables to the right and to the left. For this one should 
use the following (in fact, one of the main) assertion of the renewal theory. Let P(t, t + At) be the probability 
that the interval (t, t + At) will contain exactly one point. Then, if have density, then on the half-line the 
limit limAt^o 4rP(£,t + At) exists and tends (as t — > oo) to p — (E^)" 1 . Then the limit (as t — > oo) probability 



density of the event A(t, s), that the distance from point t to the first random point a;, > t is greater than some 
s > 0, equals 

pPr{^ =Xi- x l+ i > s) = p(l - G(s)). (1) 

That is why the first (after 0) point x\ should be taken on the random distance from with density ([1]). The 
following consecutive distances should be independent with distribution G(x). 

Alternating flows The distances in-between neighbor points are not necessarily identically distributed. The 
distributions can alternate. Take, for example, two sequences of i.i.d. random variables £1,^2, ••■ and 771 , 772, 
and put 

X2n = £l + ■•• + fn + Vl + ■■■ + Vn,X 2n -l = £l + ■•• + £n + Vl + ■■■ + Vn-1- 

Then the construction of the flow on the whole line can be done as above. 
In our case the (i.i.d.) lengths of cars and (i.i.d.) functions df alternate. 

Marked flows To each point x.- L of the point process one can assign some variable o^, taking values in some 
set S. This variable is called a mark or a spin of the point (particle, car) i, then one can call the new object 
random marked point set (flow, process). It is defined by some probability measure on the set of sequences of 
pairs (xjjCTj). Often the measure on the countable sets is given, that is the flow without marks being given, 
then at are considered to be i.i.d. random variables. In section 1.5 we shall consider marked process where the 
marks will be the velocities of the cars, moreover with sufficiently complicated correlations of the velocities with 
the coordinates. 

About Markov processes The models of time evolution of cars often use Markov processes and it is useful to 
remind some terminology. However, the mere definitions of the Markov process and its properties (for example, 
ergodicity) may differ. We give the definition we use for discrete time. 

Consider on some phase space X a system of measures (transition probabilities) P{A\x) 1 meaning the 
probability that the process at time t + 1 hits the set A C X, if at time t the process was at the state x € X. 
If all measures P(A\x) are one-point, then it is equivalent to the deterministic map T : X —> X, that is P(-\x) 
is the unit measure at the point T(x). This deterministic map defines a dynamical system. 

Note that the system of measures P(A\x) defines the mapping U of the set of all measures on X to itself 

[/> = j P{.\x)d^{x) 

The notion of invariant (with respect to U) measure is very important. Normally one studies its existence, 
uniqueness and other properties. 

Using the system of transition probabilities one can construct various sequences of random variables = x n 
with values in X, or their distributions fj, n on X, where Pr(£„ S A) = /j, n (A). The probability space is the set of 
trajectories {x n }. For example, it can be stationary Markov process £„, n G Z, where fi n are invariant measures 
on X, or non-stationary process £„, rt G Z+, with given initial distribution /iq on X. 

Markov process can be understood either as ONE of this sequences £„, or as the whole family of such 
sequences. Correspondingly, the terminology (for example the notion of ergodicity) differs. Dynamical system 
with given invariant measure /j, is called ergodic if any invariant subset has [i zero or one. Any stationary process 
can be considered as the dynamical system - time shift on the set of trajectories. Then the ergodicity can be 
defined as the ergodicity of this dynamical system. However more often, Markov process is understood more 
generally and the notion of ergodicity is defined differently. Namely, the process is called ergodic if there exists 
the only invariant measure /ionI, and for any initial measure /iq as n —> 00 we have fi n — > \l in the sense of 
weak convergence. 

Note that Markov processes normally belong to one of two classes. The first class includes all processes 
such that there exists a positive measure on X (may be infinite), with respect to which all measures P(-\x) 
are absolutely continuous. This includes all classical Markov process - finite and countable Markov chains, 
diffusion processes, etc. Such process are ergodic if there are no non-trivial invariant subsets (this property is 
called irreducibility) , and there is a unique invariant measure. In most cases convergence to it from any initial 
state follows. For countable chains an equivalent condition is positive recurrence that is finiteness (for any pair 
x, y G X) of the mean hitting time of x starting from y. 

The second class is characterized by the property that all measures P(-\x) are mutually singular. All infinite 
particle processes are in this class. The theory of such processes is essentially more complicated. 



1.2 Traffic capacity as a function of velocity 

Driver's psychology in the simplest flow Detailed modeling of the driver psychology is impossible but 
some relation is obvious. Thus, the driver i sees several cars (often only one car in front of him) in the flow and 
chooses subjectively optimal distance to the preceding car. If the velocity Vi-\(t) = rf " c ' ejt 1 ^- > of the preceding car 
changes slowly, then one can assume that the driver's reaction is faster, and the chosen distance Df depends 
of the velocity of the preceding car at the same time moment 

Dt = D+(« i _ 1 ) 

(index i says that the functions D^(v) can be different for different drivers). Let us call the flow algorithmic at 
time t, if for all i 

df(t)=Dt(vi-i(t)), 

that is all velocities are consecutively defined by the velocity of the previous cars. Of course, individual functions 
can be interesting for policemen, but we are interested with their statistical characteristics. In probabilistic 
approach the functions Df(v) become i.i.d. random variables, depending on the velocity v of the previous car 
as a parameter. The distribution of these functions cannot be deduced from mathematical or physical laws, but 
rather it depends on the individual and collective psychology of drivers and should be found experimentally, see 

M- 

Deterministic dynamics without overtaking If all cars, drivers and velocities v are identical, then many 
problems appear to be very simple. Let d be the length of cars and d + = D + (v) be the distance to the previous 
car, which the driver strictly follows. Already such dynamics allows to understand many qualitative aspects. 
Let us define the road capacity as maximally possible current through it 

J max = maxwA(i'), 

V 

where maximum is over the allowed interval of velocities, and where 

X{V) = d + D+(v) 

is the density of cars on fc-lane road for a fixed velocity v. It is clear that the flow capacity can become smaller 
when the velocity increases. This simple conclusion says only that many drivers increase the distance to the 
previous car when the velocity becomes too big. 

Random dynamics without overtaking We get the same conclusion if the velocities v are the same, but 
the functions df are random independent, and their expectations (for given v) are equal to some number d + (v). 
We see that the fact of non-trivial dependence of the road capacity from the velocity is in fact the trivial 
consequence of the driver behavior and does not need any models. However, for nicer questions stochastic 
models are necessary. Now we will introduce a sufficiently general probabilistic model with a rich spectrum 
of phases. The well-known exclusion processes appear as a degenerate particular case. Other models see in 

[E1IT01E2]. 

1.3 Random dynamics with overtaking (random grammars) 

Here we have an interesting connection with recently discovered object - random grammars, see [25j. We give 
short substantial description of one such model. 

Let at time t = all cars are situated on the left half-axis, and the movement is one-lane. We subdivide the 
lane onto intervals (cells, enumerated by Z) of fixed length and assume that there can be at most one car in 
the cell. Thus, finite sequence (group, cluster) of cars is identified with the pair (S, r), where r£Z, and S is a 
finite sequence (word) of three symbols 0,1,2 

S = SAT...S 2 Sl. 

Here corresponds to empty cell, 1 - to active (fast) driver in the cell, 2 - to a quiet driver. The length N = N(t) 
of the word and all symbols Sk(t) change with time, with the restriction that always si(t) ^ for all t > 0. At 
any time t any symbol Sk(t) has coordinate x(sk(t)). The coordinates 



x(s k (t)) =x{s!{t))-k+l 



(2) 



are uniquely denned by the coordinate x(si(t)) of the first symbol, denoted by r = r(t). 

This dynamics models the process of acceleration and slowdown of different drivers. It is defined as continuous 
time Markov chain (S(t),r(t)) on the set {(£, r)} of pairs. Jump intensities are defined as follows. Process r(t) 
models the movement of all group with common velocity v. Namely, r is increased by one with probability 
vdt for time dt, and all coordinates immediately change correspondingly to formula ([2]). Thus, S(t) describes 
replacements relative to some common motion, and is given by some random grammar, that is by the list of 
possible local substitutions (5 substitution types), that is a subword of S to another subword. Any substitutions 
of this list are independent and have different intensities (4 parameters). Here is the list: 

1. 10 — ^ 01 - fast driver moves to empty place in front of him, simultaneously creating additional free place 
after him, with probability X^dt for time dt; 

2. 120 — ► 021 - fast driver overtakes quiet driver with intensity Xjdt; 

3. 22 — > 202, 21 —> 201 - a cautious driver brakes, increasing the distance in front of him with probability 
X^dt. Note that the length S increases here (extra empty cell appears), that causes the one-cell shift of 
all cars behind. This is a non-local jump, in fact such rearrangement takes some time, but we assumed 
that this time is on a shorter scale; 

4. 200 — > 020 - a quiet driver accelerates with probability X^dt (if he decides that in front of him there is 
too much free place). 

It is necessary to say that for exact formulation of results which we only shortly describe and to get various 
qualitative types (phases) of movement, one needs to do various scalings of parameters t, N, A, We shortly 
describe only 3 phases. 

1. If A^ are small relative to other two parameters, then cars of type 2 move synchronously with common 
velocity, and fast cars have additional relative velocity. If the number of fast drivers is small, then this 
relative velocity is defined by the movement of one car among static obstacles and depends on the density 
pi of type 2 cars and density of holes po, it is approximately equal to 

Vrel = XqPO + 2A+p 2 

2. If A^~ is small with respect to other parameters (no non-local effects), and A J has the same order as Xq , Xf , 
then the difference between types disappears. We have a process close to the so called TASEP - totally 
asymmetric exclusion process. If 

\ = \ J ^2 = 
then it coincides with TASEP (about TASEP see [12]). 

3. If Aj" is small, and A^~ is large with respect to other two parameters, then the picture is different. Any 
overtaking 120 — > 021 immediately forces braking of the car 2 and, as a consequence, ALL following cars 
slow down. For the cars closer to the end of the word the slow down will be very essential, if the flow is 
sufficiently dense (few cells with zeros), as many type 2 cars will brake. 

One can complicate the introduced dynamics, for example, to avoid discretization (see the end of this section), 
introducing positive real numbers instead of zeroes - distances between consecutive cars. This demands essential 
reformulation (see section 1.6), in particular for type 3 jumps, but rough qualitative effects will be the same. 

1.4 Growth of a jam 

If the incoming traffic flow to some fixed domain equals J,- n , and outgoing flow is J out < Ji ni the number of 
cars in this domain increases in time t on 

t\Jin Jout) 

This is true, however, only if this domain is not situated on the road itself. For example, if the number of cars 
in the jam on the road increases, the answer is different.. The problem is that the domain itself can grow because 
of the arriving cars. To make this more precise one should find an appropriate model. 

Assume that the cars of the same length d move on one lane road with velocity v on the same distance d + 
between consecutive cars. During some time t the movement has been stopped by some obstacle, for example 
by red traffic light. Assume that any car stops on the distance d^ < d + from the previous car. 



Exersize 1 Prove that during time t — ¥ oo the jam (that is maximal length L(t) of the part of road, where all 
cars stand still) before the obstacle will have the length asymptotically equal to 

L(t) -t^oo tv d + d o (3) 
w d+ - d+ w 

It seems that this result depends only on the mean values and remains true even if overtaking is possible. 
This was proved in |16| for independent movement of the cars (that is when the cars do not hinder each other), 
moreover the car velocities can fluctuate but have the same values and equal v. The proof is absolutely not 
trivial. 

Local widenings and contractions of the road What occurs when a part of road with k lanes turns 
to the part with I lanes. Let this occur at the point x = 0. 

Case k < I. Assume that maximally allowed velocity equals v max and the drivers are assumed to be 
disciplined. All cars move along fc-lane road with velocity v < v max . And it is impossible to move faster due to 
the fundamental relation 

d + D+{v) = p- 1 

between density p of cars and their velocity. Then along /-lane road the cars could conserve the density and 
move with the same speed, but p can change so that the cars will move with some greater speed v\. The time 
gain is 

L L 

V Vi 

Case k > I. Here three different situations are possible. 

Free flow If the flow is very sparse, then the cars arrive at the point alone and will not notice the decrease 
in the number of lanes. 

Growing jam Denote Jf. the incoming current and let Ji. ma x be maximally possible current along /-lane 
road. If Jk > Ji,max, then the jam will grow, and the number of cars in this jam will grow as t(Jk — Ji,max), 
more exactly as in formula ([3]). 

Delay For the case Jk < Ji.max practical observations say that at point jams of random lengths can 
appear, which however will not grow too much. There are no corresponding stochastic models. To create such 
model, it is useful to have a collection of models for shorter time phenomenon, Namely, when the standing cars 
at the traffic lights start their movement. We shall describe some models of this type now. 

1.5 Speeding up 

In [5] the cars are points 

... < x l (t) < Xi--y{t) < ... 

on the line. At initial time t = the cars stand still and are described by Poisson point field with density 
p < 1. The cars can have two velocities: and 1, and overtaking is prohibited. Any incumbent car starts its 
movement with speed 1, independently of the others, in exponential time with mean 1. It can occur that the 
car with number i will reach the car i — 1 when it did not start its movement. Then the car i stops and resumes 
its movement in exponential time after the moment when the car i — 1 starts movement. This rule acts at any 
time. The model describes roughly how clusters of cars start to move. 

The main result is that with probability 1 each car will stop finite number of times (if p < 1). Let U be 
the minimal time moment after which the car i does not stop anymore. Then for any i and k and any time 
moments ti,ti_i, ...,ti_fe the random variables 

^i-i(^-i) - Xi(U), Xi-2(ti- 2 ) - a;i_i(ti_i), ...,Xi-k(ti-k) - Xi-k+i(ti-k+i) 

will be independent and exponentially distributed. Otherwise speaking after leaving the jam, the cars will form 
Poisson configuration with the same density p. 

Assume now that at time all cars are situated on the left half-axis. The distribution is again Poisson with 
density p. Each point moves with velocity v > 0, if the distance to the next car from the right is not less than 
some d e ff > 0, and stands still otherwise. Here it is evident that any particle will not stop starting from some 
time moment. But in this model one can get more. Consider the following random variables: random time t£ 

(2) 

when k-ih point starts movement, random time T k when this point does not stop anymore, the distance Xk 

(2) 

between points k and 1, starting from time moment tv . 



Problem 1 Find asymptotics of distributions of these random variables when k — > oo. 

The connection with delay problem is evident. Let we have two lanes and on any lane the flow density be 
p; the joint flow thus have density 2p. The cars from the first lane now want to squeeze into the second. The 
squeezing algorithms can be quite different. For example, any car can squeeze independently of the others, if 
its distances (along x-axis) to the cars from the second lane are not less than some d + . 

1.6 Short and long range order with time dependent velocities 

Here the cars are represented by points Xi. To the car i a stationary random process Wi(t) is assigned, which 
defines time evolution of its velocity on empty road (that is when there are no obstacles in front of this car). 
This process implicitly defines the activity of the driver at given time t. The processes u>i(t) are mutually 
independent and are defined only by the psychology of the individual driver. Assume that there are constants 
< C\ < C2 < 00 such that for all t, i 

d < Wi {t) < C 2 . 

The flow is given by initial coordinates Xi(0) of the cars, and their movement is defined as 

Xi(t) = Xi(0) + / Vi(s)ds, 
Jo 

where Vi(t) is the velocity of i-th car, defined below. Moreover, we assume that the initial coordinates are such 
that the distances df(0) are independent and, for example, are exponentially distributed with parameter p(0). 
We shall say that the car i has an obstacle in front of itself at time t, if 

Xi(t - 0) = Xi-i(t) 

The process will be completely defined if for any t\, t n , i\, i n we could define finite-dimensional distribution 
of the vectors 

(Vi^ti), ...,v in (t n )) 
where some ik can coincide. The following rules do this job: 

1 . (rule for free road) if none of the cars i 1 , . . . , ij~ for k < n does not have in front of itself an obstacle , then the 
distribution of the vector v il (t\), Vi k (tk) coincides with the distribution of the vector Wi 1 (t\), Wi k (tk) 
and is independent of the distribution of the vector Vi k+1 (tk+i), ...,Vi n (t n ); 

2. (rule for an obstacle) if the car i has obstacle in front of it at time t, then Vi(t) — Uj_i(t); 

3. (rule for overtaking) if the car i has obstacle in front of it at time t, then it overtakes this obstacle with 
some intensity A, during the whole (random) time interval while Wi(t) > Vi-\(t). The meaning of this 
condition is that the driver overtakes if his activity is higher. 

Even for such simplest behavior of drivers and corresponding definition of the flow there are many problems. 
We formulate some of them. 

Let us call free phase the case when the intensity of overtaking equals infinity. Then for any cars their 
velocities are independent and the covariance is zero, that is 

Covij(t) = Evi(t)vj(t) - Evi(t)Evj(t) = 

The complete order phase is when A = 0. Then there is no overtaking at all and the cars move as a queue. 

Problem 2 Assume that all Wi{t) have the same distribution. Is it true that for small A the movement will 
resemble complete order phase, and the large A will resemble free phase ? Does exist some third phase (cluster 
phase) for intermediate values of A ? 

Possibly to get more interesting qualitative picture one has to define the same process with the lengths 
di,df, and with additional indices corresponding to several lanes. Also the behavior of the driver can depend 
not only on the car in front but also of the car just behind. 



1.7 About relation of stochastic approach with kinetic equations 



For any particle xi on the real line there is a mark - velocity v\ > 0. Then the flow of cars is defined by countable 
subset {(xi,V{)} of the phase space R x R + . And the random flow, at a given time, is a probability measure /j, 
on the set of such configurations. Let n(A) = n fl (A) be the mean number of particles in bounded set A of the 
phase space. If there exists a function f(x, v) such that for any A 

n(A) = / f(x,v)dxdv, 

J A 

then / is called one-particle (correlation) function of the measure \x. 

Although kinetic equations approach to transportation problems is known long ago [33] and was studied 
a lot in physical papers [34, 35] we think that rigorous mathematical deduction of kinetic equations (that is 
equations for f(t] x, v)) from dynamics of individual cars stays open. Here we give only short comments. 

Let at time t = be given a distribution /i with one-point function f(0, x, v). If particles move freely, that 
is each particle moves with fixed velocity, different for different cars, then 

f(t + S; x, v) = / (t; x — vS, v). 

Subtracting /(t; x, v) from both parts of this equality, dividing by 5 and taking limit 8 — > 0, we get 

m +V dx-=°- (4) 
This is the simplest kinetic equation without collisions. Its solution is of course 

f(t;x,v) = f(0;x-vt,v). (5) 

In general kinetic equations should contain collisional term. What is collision in traffic ? It could be 
close approaching of cars, that is the reason for the drivers to change their behavior, for example to perform 
overtaking. Also one can consider formation of a new compound particle - cluster of several cars with short 
lifetime, which decays shortly on smaller clusters or separate cars. Moreover, there can be clusters on various 
scales - such situation had never been formalized even in mathematical physics. 

Another possible approach is to consider movement of fast cars in the media of slow cars. Let there exist two 
types of cars with velocities v and V\ > v correspondingly. Slow cars move independently and their one-particle 
function looks like ([5]). In the coordinate system moving with speed v, slow cars will become standing obstacles. 
Rules of overtaking define velocity v\{x) of fast cars, in the new coordinate system, in the vicinity of separate 
obstacle. This can give, under appropriate scaling, some kinetic picture. But as far as we know there were no 
efforts to write down rigorously even such simple models. 

There are however papers where Burgers equation is deduced from stochastic particle picture. However, the 
corresponding particle process is far from one-dimensional car flow, see reviews on exclusion processes |31[ 132] . 



2 Calculation of mean velocity on the road 

We consider here simplest problems concerning diminishing of the road capacity because of random fixed 
obstacles (accidents, repair work) and because of slow cars. Our goal is to show (solving completely some model 
problems) that it is possible to get simple formulas which allow to understand main reasons why the mean speed 
can decrease. Our main assumptions is the homogeneity of the road, exits and gates and other parameters. 

2.1 Road as one-dimensional queuing network 

The following model is taken from [8], p. 117. Let we have infinite road and two types of cars, given by points 
of the real line, the points move in the right direction. First type (faster) cars move with constant speed V\, 
cars of the second type (slow cars) have constant speed V\ > v 2 . 

Assume that faster cars initially (at time t — 0) have Poisson configuration on R with density Ai. Slow cars 
initially are at the points 

x = < X\ < ... < x n < 

where the distances Xk — Xk~i are identically distributed with mean A^ 1 (not necessarily exponentially). Slow 
cars move completely independently - other cars do not influence them. Faster cars «interact» with any car 
when their coordinates coincide. Namely, they can overtake slow cars. When a faster car catches up with a 



slow car, that is their coordinates coincide, then it follows slow car for some time, that is moves with the speed 
V2- In exponential time with parameter fi the fast car overtakes the slow car and moves again with the speed 
v\. If the slow car catches up the group of fast cars still following some slow car, then the overtaking occurs as 
in a queue, in order how the fast cars follow the slow car. 

Without loss of generality one can put Vi = 0, and the speed of fast cars put equal to v — v\ — v-i- That 
is why each slow car can be considered as a server, where the clients (fast cars) arrive and are waiting for the 
service (overtaking), they are served with intensity /i. 

Now the problem can be reduced to the linear queuing system which we shall describe. There is infinite 
sequence 

5*0 - ► ••• — >• Sk — > Sk+i — > ■■■ 

of nodes (servers) of two types. Each server Sk is a M/M/l type server with FIFO (first-in- first-out) discipline, 
that is in the natural order. These servers correspond to slow cars, and clients correspond to fast cars. For 
example So corresponds to the extreme left slow car. Second letter M means exponential service rate. First 
letter M means Poisson arrival stream. Thus the clients to So arrive as a stationary Poisson stream with 
intensity X\V. From elementary queuing theory it has been known, firstly, that if X\v < /z, then the stationary 
regime will be established with the probabilities 

P n = (l-r)r",r=^. 

M 

that the length of the queue equals n. Secondly, (Burke theorem), that in the stationary regime the output flow 
from M/M/l will be Poisson with intensity equal to intensity of the input flow, which is X\v in our case. 

After the first node, the output flow, with random but the same for all clients time shift Xl ~ x ° , arrives to 
node Si, where also stationary regime will be established. 

Let us find mean velocity of fast cars on the interval (x ,xn), N — > oo. We shall assume that the stationary 
regime has already been installed. The time along this path is the sum from N overtaking and N paths between 
slow cars. Mean overtaking time is 



g ( i_ rK (n±i) = i 



n=0 

and mean time to catch up the next slow car is 
Thus, the mean speed of fast cars is 



(1 — r)ji jj, — Xiv ' 



1 

X2V 



Xn 1 



~" tea " ( M -A lW )- 1 + (A 2W )-i- 
In the next section we will consider more complicated situation with more general distributions. 



2.2 Mean speed slow down due to repair works 

On a long road the cars move with constant speed v, but encounter obstacles. The obstacles have small size 
comparative to distances between them, and we assume the obstacles are just points. They can appear on 
arbitrary part of the road (a;, x + dx) C R during time interval (t, t + dt) C R with probability Xdxdt. More 
exactly, the pairs (coordinate and time moment of obstacle appearance) (xj ,tj) £ Rx R+ form Poisson random 
filed II on R x i?+with intensity A. Another equivalent definition is that for any time interval I C R there is 
Poisson arrival stream with intensity A|ij, moreover the arriving obstacle choose the point uniformly on the 
interval /. 

Assume that j— th obstacle stays on the road for some random time tj, after that it is cleared off from 
the road. Random variables tj are assumed to be i.i.d. with distribution function Q(t) and independent from 
Poisson random field II. Assume finiteness of mq = Etj and = Etj. 

Further on we consider two cases. In the first case the bypass is prohibited and the car stands until the 
obstacle will be deleted, after this the car immediately moves with speed v. In the second case, assume the new 
car arrives to an obstacle and sees n cars, which have not yet bypassed this obstacle. The new car is allowed 
to bypass the obstacle and these n cars simultaneously. This takes takes some random time which does not 
depend on n. Denote rj m ^ the random bypass time of m-th obstacle by i-th car. Random variables rj m ^ are 



assumed to be i.i.d. with distribution function F(u). These assumptions are natural for small density of cars, 
with small number of cars accumulating after the obstacle. Below we consider the case of greater load. 

First of all we shall find the mean speed of the car. Under above assumptions the cars do not impede each 
other, that is why it is sufficient to consider one car problem. Denote T(x) random time the car passes distance 
x. We will find the limit of the ratio as x — >• oo. 

Let b = AmQ, and £ be a random variable with density 

h(t) = m Q 1 (l-Q(t)). (6) 

Note that 

oo oo oo (2) 

EC = — [t(l - Q{t))dt = — f(l - Q(t))d (?) = J- / t 2 dQ(t) = ^_ 
m,Q J rriQ J \ 2 J 2m Q J 2m Q 

oo o 

(2) 

where ttiq is the second moment of Q(t). 

Put a — min(77,£), where equality is in distribution. Moreover, the random variables n, £ are considered 
independent and rj has distribution function F(u) . Put 

a = Ea. 

Theorem 1 With probability 1 as x — s> 00 

x v 
T(x) 1 + abv' 

Proof. With no loss of generality assume that the car enters the road at the point x 
Tg(x) be the idle time of the car. Then obviously T(x) ~ Tq(x) — v~ 1 x and 

xx 1 
T{x) = T{x)-T {x)+T (x) = v-^+x-^T (x) 

Thus it is sufficient to find the limit of the ratio as x — > 00. We want to show that 

?r(x) 

t q( x ) = ^2 a " 

i=l 

where are i.i.d. random variables distributed as a, 7r(x) is a random variable with Poisson distribution with 
parameter bx. Random variables oti and n(x) are assumed to be independent. The meaning of this formula, 
that the car on the distance x will meet tt(x) obstacles and will loose random time on i-th obstacle. 

From ((8J and strong law of large numbers it follows easily that — > ab a.e. as x — > 00. 

Let us prove ([5]). Introduce marked Poisson point field IL on R x R + with configuration (xj,tj,Tj), that is 
Tj is the mark at the point (Xj,tj). The following assertion can be found in [S]: 

Lemma 1 Marked point field Hi is equivalent (in probability) to the Poisson field on R x R 2 : with intensity 
XdxdtdQ(t). 

On picture 1 the obstacles are presented as horizontal segments. The coordinates of the initial point define 
place and time of its appearance (pair (xj,tj)). The length of the segment is the time of obstacle stay on the 
road (mark Tj). 



pic.1 




(7) 

at time t = 0. Let 



(8) 



Consider any straight line c\t + ci and consider its intersection points with horizontal segments. Denote 
{xi} the space coordinates of these points as it is shown on picture 1. The next lemma is proved [7]. 



Lemma 2 Configuration {xi} is the Poisson process with intensity b = Xttiq. 



Picture 2 shows the trajectory of the car, which starts from point x — at time t = 0. Denote Xi the space 
coordinates of the obstacles appearing during movement of this car, ti are the moments of their appearance, 
Si - the moments when the car meets the obstacles, Ui - time moments when the car gets rid of them either 
bypassing or because the obstacle disappear; on — Ui — Si - car delay i-th obstacle. 




pic. 2 



From lemma [2] and from space-time homogeneity of Poisson field II it follows that the points Xi make up 
Poisson process of intensity b. 

The residual lifetime of an obstacle is its lifetime after the car catches it up. In other words, it is equal to 
the car delay on this obstacle. 

Lemma 3 The residual lifetime has the distribution with density ft-(s), where h(s) is defined by (OJ). 

In fact, from the properties of Poisson point field it follows that the conditional distribution of residual 
lifetime, under the condition that the full lifetime equals t, coincides with the uniform distribution on the 
interval [0,t\. By lemma [5] the probability of obstacle appearance on the interval dx equals XniQdx + o(dx), and 
the probability of appearance of an obstacle with full lifetime t in the interval dx is equal to XtdQ(t)dx + o(dx) , 
that follows from lemma [TJ As 

XtdQ(t)dx + o(dx) _ tdQ(t) 
XmQdx + o{dx) itiq 

is the conditional probability of appearance of an obstacle with lifetime t, then the density of residual lifetime 
is 



DC 

/ 



tdQ{t) ds 



t 



771q 1 (1 — Q(s))ds = h(s)ds 



Lemma is proved. 

If bypassing is possible, the car will lose time equal to minimum of bypass time and residual lifetime of the 
obstacle, that is Ui = min(?7, £). Theorem is proved. 

Let us discuss the result. We have seen already the meaning of the constant a, the constant b characterizes 
stationary density of obstacles in the space. 

This result is rather exact for small density of cars, as there will not more than one car at the obstacle. For 
large density of cars the bypass time will increase dependent on the mean length of the queue at this obstacle. 



2.3 Mean speed slow down due to slow cars 

The road here is the real line R. Flows are not too dense, thus the length of cars does not play role, and at 
a given moment the position of the car is given by the point Xi(t) £ R, where i is the index enumerating the 
cars. Each car has predefined route: place and time of entrance a^.m, t^m, and the exit coordinate Xi jOU t- But 
the exit time ii iOM t depends on the road capacity. We define mean velocity of the car i as 

Xi^out Xi^i n 

ti : out ti^iri 

There are two types of cars - fast and slow, any car moves with constant speed in the right direction. Fast cars 
have speed vi, slow cars have speed t>2, where v\ > i>2 > 0. Put v = v\ — V2- Slow cars move along the road 
without stops, and fast cars - until they reach a slow car. After this the fast car i moves with slow car j for 



some random time Tjj, and them bypasses it immediately getting speed v\. Main assumption is that random 
variables r^.j are i.i.d. with distribution function F(s). 

This distribution function may be found statistically by two ways - by direct sampling and by estimating 
the density of obstacles, for example for one-lane road it can be the density of the opposite flow. 

Arrivals of slow cars is given by the same Poisson point field IT with intensity A, which was defined in the 
previous section. We will need more notation. To any slow car j we assign random distance pj which it should 
pass, after that it will leave the road. Random variables pi are assumed i.i.d. with distribution function G(r). 
Moreover pj does not depend on the random field IT. Assume existence of two first moments itiq = Epi, 

™( 2 ) — J7 n 2 

m G — h,p x . 

Slow car does not meet obstacles and passes all its route with speed v%. Fast cars can be stopped by slow 
cars. Consider two cases. In the first case the overtaking is prohibited and fast car follows the slow car until the 
latter will leave the road, after that the fast car immediately gets its speed v\. In the second case the overtaking 
is possible. More exactly, when i-th fast car catches up j-th slow car or a group of fast cars following j-th slow 
car, it takes random time Tjj to bypass j-th slow car or all group of these cars. This time does not depend on 
the size of the group. Tij are assumed to be i.i.d. with distribution function F(u). 

Put d = \ma (v^ 1 — v i 1 )- Introduce random variable (3 with distribution density g(x) = mg'(l — G(x)) 
and put 7 = min^T^i, /3), where the equality is in distribution, and random variables t\,i,P are assumed 
independent. Note that 

Put c = £7. 

Theorem 2 With probability 1 as x — > 00 

x _ 1 + dc 

7f7~T v i = -r—j —vi- (9) 

±{x) 1 + dcviv 2 

Proof. Let us note that this case can be reduced to the case v 2 = 0. Introduce new coordinate system moving 
with the speed v 2 relative to the initial system. Find the mean speed of fast cars in the new coordinate system 



EP 



(2) 
ra G 

2m G 



with the formula ([7]), substituting v — v± — 112, b — " G , a — ^ 



V2 ' V2 
I Vl - V 2 



(v x - vz)- 1 + ^ l + dcvxv- 1 

v 2 

Then the mean speed of a fast car relative to the initial coordinate system is 

vi — vi 1 + dc 

vi — rr + V2 = rr^i- 

1 + dcv\v 2 1 + dcv\v 2 



3 Analysis of complex transport system 

Graph is a natural tool to describe transport networks (for example, city streets), where the set V of vertices 
represent crossroads (nodes or service stations), and the set L = {(«, j)} of edges represent segments of streets 
without crosses. Let N be the number of vertices. We assume that there not more than one edge between two 
vertices. 

Mostly the following two classes of networks were studied. First class includes (by the name of the authors 
and in the order of increasing generality) Jackson networks, BCMP and DB networks (see [S], [2S])- Therein 
the client (communication, car, job), after being served in some vertex, chooses randomly the next vertex. The 
second class are the Kelly networks (see [S]), where each client has apriori fixed route. These classes have much 
in common - close results and techniques. For example, they may have very useful multiplicativity property, 
that is the stationary distribution has the so called product form. Here we consider only the first class of 
networks. 



3.1 Closed networks 

If the cars do not come from without and do not leave outside, the network is called closed. Then the number 
of cars in the network is conserved and we denote it by M. The movement of a car is defined as follows. The 
car is waiting for some time on the cross i and is directed afterwards to some vertex j. The choice of j is defined 



by some stochastic matrix, we shall call it the routing matrix, P = {pij}i,j=x,...,N) where pij is the probability 
that from vertex i the car will go (after service time) to vertex j (for example, to the right, left or straight on), 
that is along the street 

Stochastic matrix P defines finite discrete time Markov chain with state space V — {1, N}. We assume 
this Markov chain irreducible. Then the system of linear equations 

N 

PP = p, p= {p!,...,p N ) ^PiPij = Pj, j = l,-;N (10) 

2 = 1 

has a unique solution (up to a common factor). Normed solution gives its stationary distribution 

Pi 



Li=l Pi 



1,...,N 



For each vertex i 6 V we define a function p,i(rii) of the number rii of cars in i-th vertex, where Pi{0) — and 
Pi{ni) > for rii > 0. This function characterizes capacity of this node and defines the intensity of output flow 
from this vertex. Namely, the probability that from node i, during time period t, exactly one car will leave the 
vertex is pi{rii)dt + o{dt), under the condition that in the node there are rii cars at time t. Using queuing theory 
terminology we call pi(rii) the service intensity in vertex i. 

The order in which the arriving clients are served is defined by service discipline (protocol). The simplest 
possibility is FIFO - first come first served. If there are n, cars in the vertex i, then the first car in the queue 
is served with intensity //j(rij). 

More general service discipline is the resource sharing discipline, the resource here is the capacity of the 
crossroad. This means that resource is divided in some proportion between cars in this crossroad. We assume 
that that fc-th car in i- th node is served with intensity /ii,fc(n,') < Pi{rii) so that 



fc=i 



P-i,k( n i) = Vi( n i)- 



For example, common resource can be divided equally between each car in the queue 

, > Pi{n r ) 

Pi,k{ n i) = 

Tli 

If Mi.i( n i) — Pi( n i)i then we get FIFO discipline. Thus, the intensities pi.k{ni) completely define service 
discipline in the nodes. 

Dynamics in the network is described by N— dimensional continuous time Markov chain £(t) — (£i(t), i = 
l,...,N), where £i(t) is the number of cars-at vertex i at time t. Thus, the state space Sm is the set of all 
vectors with nonnegative integer coordinates n — (ni, Ujv), such that n% + ... + tin = M. 

Let ej be a base vector, where «-th coordinate is 1, and the rest coordinates are 0. From state n Markov 
chain £(t) can jump to the state T^n — n — + ej,i ^ j, with intensity 

a{n,T itj n) = Pi(ni)pi t j, (11) 

under the condition that rii ^ 0. The jump n — » Tijfi corresponds to the fact that after leaving cross i the car 
arrives to cross j. 

Note that Markov chain £(t) is uniquely defined by the routing matrix P and service intensities (p,i(ni),i = 
1 V). 

Let p = (px, pn) be the solution of equations (|10[) . which we consider as a formal equation for intensities 
Pi of input flows to the nodes (in the stationary regime these are equal to output intensities). Solving these 
equations, we find pi. Then the stationary distribution v(ni, ...,npf) of Markov chain £(t) looks like 

1 N 

u(m,...,n N ) = - 1] m ik r~v ( 12 ) 

/jv,m p l {l)p l {2)...pi{ni) 



where the normalizing factor (the canonical partition function) 

N 

' J N,M 



E II w (l) w (2)... w (n<)' 



ni + ...+njv=Af i- 



This can be verified by direct substitution of (|12[) to Kolmogorov equations for stationary probabilities, see for 
example [S]. 



3.2 Open networks 



Consider network consisting of N nodes. Contrary to close network, total number of cars is not fixed. Assume 
that to node i arrives (from outside) Poisson flow with intensity Aj, i G {1, N}. 
The routing matrix P — {pij}i,j=i n is assumed to be decomposable and 

N N 

Vi : J2 p v - 1; 3i ° : J2p^ < L ( 13 ) 

Similar to closed network, pij is the probability that a car from node i goes to node J. The difference is that 
there is also probability 

N 

Pio = 1- 5^P»i- 

i=i 

for a car at node i to leave the network. 

As for closed network, let pi{iii) be the service intensity at node i. Then the car at node i leaves the network 
with intensity pi(rii)pifi. 

The dynamics of an open network is described with TV-dimensional continuous time random vectors rj(t) — 
(r)i(t), i = 1, N), where r]i(t) is the number of cars at node i at time t. Random process rj(t) is a continuous 
time Markov chain with state space S, where S is the set of iV-dimensional vectors with non-negative integer 
components n = (m, njv)- From the state n Markov chain can jump to one of the states T^ n = n—ei+ej, 
Tififi = n — ei, Tifi = n + with intensities 

a(n,Tijn) = ^{n^p^j, 

a(n,T lfi n) = ^("-i)Pi,o, (14) 
a{n,Tin) = A 4 , 

under condition that Tijh, Tj^n, Tifi G S. 

Thus, Markov chain r)(t) is uniquely defined by the triplet (A,/i,P), where A = (Ai, ...,Ajv) is the vector of 
intensities of input flows, /i = (/ij(nj), i = 1, N) - the vector of service intensities in the nodes and P is the 
routing matrix. 

Consider formal equation for input intensities to the nodes (in the stationary regime they are equal to output 
intensities) 

N 

p = \ + pP <j=^ pi = Xj + PkPH, Vi (15) 

fe=l 

Under the condition (|13[) and the irreducibility of P, this equation has a unique solution which can be written 
as follows 

oo 

p = A + AP". 

n=l 

Consider now the case when the service intensities p,i{rii) = pi do not depend on the number of cars in 
nodes. Define the loads in the nodes by the formula 

Ti = — , i = 1, N. 
Pi 

The following theorem can be found, for example in [3], it is called Gordon-Newell. 

Theorem 3 Markov chain rj(t) is ergodic if and only if for all i = 1, ...,N we have Ti < 1. Then the stationary 
distribution is given by 

N 

a{n^...,n N ) = - n)r"*. 

i=l 

It easily follows from this theorem that the mean queue lengths in the stationary regime are- 



III; 



1 



If in some nodes i±, the load is strictly larger than 1, then Markov chain?7(t) is transient. Then the mean 
queue lengths at nodes i\, ...,ik tend to infinity as t oo. A detailed analysis of open networks see in [TS]. In 
particular, it is shown that in the nodes with the load exceeding 1, mean queue lengths grow linearly with time. 
There are also explicit formulas for the mean speed of growth. 



3.3 Algorithm to find the critical load in closed networks 

This section is based on the paper |18] . Consider a sequence of closed networks J/v, N — 1,2, .... Network Jn 
consists of N nodes and M = M(N) cars. Service intensities at the nodes of network Jn do not depend on the 
queue length : fii ff{ni) = Hi,N- Let Pn = {Pi,j,N} be the routing matrix in the iV-th network; Pn is assumed 
to be irreducible. 

Let pn = (pi,N, ••• ) /9iV,iv) be the vector with positive components satisfying equation 

Pn = PnPn- (16) 

Relative loads in the nodes are defined as 

Ti,N = C N Pi t NTi,N, 

where t^n = (\n an£ l — m&Xi=i,...,Ar Pi,NT~i.N- It is evident that r^N & [0, 1]. 

Correspondingly to (|12|l the stationary distribution of the number of cars £i,jv,M m the nodes of J/v is equal 

to 

1 N 

Pn,m(&,n,M =n l ,i = 1, N) = — TT r i!jv 

' 2—1 

where the normalizing factor (canonical partition function) 

z N , M = n<v- (17) 

n\ + ...+nM—M i=l 

Many important network characteristics can be expressed in terms of this partition function. 
Exersize 2 Show that the mean number of cars in the i-the node in the stationary regime is 

Pt r i,N dZ N<M 

m>i,N,M = &t,i.N,M — — ^ 

AN,M Ori^N 

We are interested in cases when N, M are large enough, more exactly N, M — !• oo and so that j£ — >• A = const, 
that is the number of cars per one node is constant. We shall see that A defines existence or non-existence of 
jams in the network. 

To formulate results as general as possible, we shall assume weak convergence of relative loads ■ More 
exactly, define sample measure on the interval [0, 1] 

where A is arbitrary Borel subset of [0, 1]. Assume that as N — > oo the measures In weakly converge to some 
probability measure I on [0, 1]. 

Remark 1 It could be interesting to find concrete sequences of growing networks for which the limiting measure 
I can be found explicitly. Some examples where measure I is one-point one can find in the bibliography to the 
paper ]18]j . 

In terms of limiting measure I we shall find critical value X cr of the density, so that for A < X cr mean lengths 
of the queues are uniformly bounded. If A > A cr , then at the node with maximal value of load the mean length 
of the queue tends to infinity, that means the existence of a jam at this node. Put 

where z G C \ [1, +oo). The function h{z) is strictly increasing on [0, 1). Put 

A cr = lim h(z). 
z— >i— 

We shall assume that A cr > 0. 



Theorem 4 • If A < X cr , then mean lengths of queues are uniformly bounded: there exists constant B, 
such that mi t N < B, uniformly in N > 1 and 1 < i < N . 

• If A > X cr and i(N) satisfies condition rj(jv),jv = 1; then m^jv).^ — > oo, as N — > oo, t/iat is jams wiiZ oe 
in the nodes where the load is maximal. 



ForzeC \ [1, +oo) put 



1 N 

S N {z) = -A(l + e N ) In z - — ln (l - zr i>N ), (19) 



AT 



S(z) = -Xlnz- / ln(l - zr)dl(r) 



where A(l + e w ) = M. 

Introduce the generating function (the grand partition function) 

JV 



e N (z) = ]r z M ^ M = n 1*1 < 1 



M=0 

By Cauchy formula and (ITU|) we have the following expression for the partition function (IT71) : 

27TI J 27TI J Z 

7 7 

where 7 = {z € C : |z| = cr< 1}. For the means, accordingly to (fl"8|) . we have 

mi.jv = / 1 exp(NS N (z))dz. (21) 

27rzZjv 7 1 - zr hN 
7 

One can show that for the stationary distribution of queue lengths the following formula holds 

1 f _ n 

Pn,m(£i,n,m = ni, ...,£k,n,m = n K ) = - .„ / z^ 1 TT(1 - zr iiN )(zr iiN ) ni exp(NS N (z))dz. (22) 

7 i=1 

In the proof of theorems of this section, essential role is played by the steepest descent method (see, |27J), 
more exactly, its generalization concerning the case when the function in the exponent depend on N. From 
equation 

^1 = (23) 
oz 

we find steepest descent points. Let zq,n be the root of this equation lying in (0, 1). 

Exersize 3 Show that all roots of equation i23\) are real and positive. Always there exists unique root, which 
belongs to the interval (0, 1). 

Let Zo-be the root of equation 

, . . A dS(z) , . 

lying in the interval (0, 1). 

Exersize 4 • Prove that for all A there exists limit limjv-»oo Zq.jv = z a = z o(X) > 0. 

• If A < A cr , then zq(X) is the root of equation {2$ ; Zq{X) is strictly increasing in A, Zq(X) G (0,1), 
lim A ^A cr _ z (X) = 1. 

• // A > X cr , then Zq = 1 . 

In the next theorem we find the asymptotics of the partition function and limiting distribution for the 
sequence of closed networks Jn- 



Theorem 5 Let A < X cr , then 

• As N — > oo the partition function Zn and the free energy Fn — In Zn have the following asymptotics 

exp(NS N (z , N )) 1 

• If for i = 1, ...,1^ £/iere exist i/ie limits n = limAr_ ! . 0O r^jv, i/ien 

lim mijv = , 

jV->oo 1 — Zofi 

K 

111 _P;v Af (£l MAT = Tit , .... /V Af = tik) = 
TV— >oo 



lim Pn,m(€i,N,M = ni, -,£k,n,m = n K ) = TT(1 - z r i ){z r i ) r 



Thus in the limit we get an open network consisting of independent queues. 

Proof of theorems [4] and [5] We give more general result from which the theorems [4] and [5] follow. Let 
Ud(v) = {z G C : \z — v\ < d}. Consider the contour 7 = {z G C : \z\ = zo(X)}. 

Theorem 6 Let A < A cr and f(0,z), 8 <E Q, be a family of functions holomorphic in the ring {z G C : 
zq(X) — So < \z\ < zo(A) + <5o} for some 5a > 0, uniformly bounded in this ring and such that for sufficiently 
small e > there exists such 5 U > and nonzero constant f u , such that \f(8,z)/f u — 1| < e for z G U26„(zq), 
6 G 8. 

Then, for N sufficiently large, uniformly in 8 G 

1 f tta \ iatg < \\a f u exp(NS N (z , N )) n . r > 

— / f(6,z)exp(NS N {z))dz = — 1 + Cw , 

2tti y J2-kNS"(z ) 

7 

where \C,n\ < 25e. 

Proof of this theorem is based on the saddle-point method, see [27]. Difference from the standard situation 
is that the function in the exponent depends on N . Detailed proof one can find in the original paper [18J. 
Proof of theorem [5j Using theorem |6] we prove the first part of the theorem [5j Using ([20]) we have 

1 f exp(7V Sn(z)) 

<6n,m — - — ; / az, 

Itxi J z 

7 

where 7 = {z G C : \z\ = 20(A)}. Putting f(0, z) = z~ x , f u — z^ 1 and using theorem [51 we get that for any 
sufficiently small e > for sufficiently large N 

exp(NS N (z 0tN )) 

Zn = 7S=^77^^( 1 + l&H < 25e - ( 25 > 

zoV 2 ^^ (z ) 

The second part of theorem [5] can be proved similarly by using (|2Tj) for the mean queue length and (|22|) for 
joint distribution of the queue lengths. 



Exersize 5 Prove third part of the theorem [3J using theorem [3| and formulas \21\i 

Proof of theorem 21 To prove the first part of theorem HI consider the family of functions 

f(8,z) = - + -?—, 969 = 10,1], A>0, /„ = -. 

z 1 — za zq 

Fix small e > and choose a u = | , A = ( 1 16 ^°- )e ■ By theorem [6] we have for sufficiently large JV and all 8 G 

1 f f A 8 \ . ^ , ,, , Aexp(NS N (z Q N )) . , |A . „_ 
M >exp(iV^(z))dz= 1 J!l " / ■(! + Cat), ICjvI < 25e. 



2mJ \z 1-zOJ " z y / 2TrNS"(zo) 



Dividing by Zjq and applying (|25[) to the right-hand part of the resulting equality, we get for sufficiently large 
N 

A+^^- f -?—eMNS N (z))dz = A(l + ( N ), \( N \ < 30e. 
Zn Zm J 1 — za 

7 

From last equality and formula (|2ip we have uniform boundedness of to^at. 

Prove the second part of theorem [4] For this we need the following monotonicity property: for any Mi > 
Mi > and any N > 1 the inequality mi,M 2 ,N > Tni,M t ,N holds. 

As zo(X) is strictly increasing in A, zo(A) G (0, 1) and limA->A cr - Zo(A) = 1, then the function 

2o (A) 



1-*) (A) 

is monotone increasing and tends to oo, when A /*■ X cr . That is why for any m > there exists such A' = 
A'(m) < A cr , that 

1 - zo(A') 

Without loss of generality we can assume that i(N) = 1 and r± t N = 1- If we put M'(N) = [A' AT], then by 
theorem [5] 

It follows that for sufficiently large AT 

1-30 (A') 

But M/N -> A > A cr > A', that is why for sufficiently large N we have M(N) > M'(N). By monotonicity 
Tii,Ar = mi,M(N),N > toi,m',at > for sufficiently large AT. This proves that mi,jv — > oo. 

Technical generalizations and mathematical problems We assumed above instantaneous displacement 
between crosses. I.e. we did not take into account time spent along the streets. This can be easily overpassed 
by introducing more extended graph. Namely introduce additional vertices uy, corresponding to the street 
between crosses i and j, and mean duration Tij=p^ 1 of movement along the streets. In the queuing terms 
this means that the streets are considered as new service nodes with infinite number of servers and exponential 
service time with mean Tij=p~ J 1 . 

Note that the results of section 3.1 can be generalized to the case, when the network contains nodes with 
infinite number of servers. Let, for example, the network contain one such node (i = 0) and po,N{n) = riv^ 
is the service intensity in this node. Let p^ — (fo,N, ■■■,Pn,n) be the solution of equation (|16l) . Then relative 
loads can be defined by formula 

H0,N pi.N 

n,N = ) 

P0,N fJ*i,N 

so that ro,jv = 1- According to (fT2"| . the stationary distribution of the queue lengths is 

1 1 N 

Pn,m(&,n,m = n h i = 1,...,N) = ^ M wII r Dv' 

Zn.m {M -2Z i= iiHV- i= i 

where 



z n ,m= — ^Il r 



and the grand partition function is 



N 1 

AJL 1 - zn N 



N 

"i,N> 



Put 



q i N = , w = zp N , p N = max r i N . 

PN Ki<N 



Then 

N 1 

Under assumption that pnN — > a > as A^ — > oc we can find critical value of density A with the formula 

l 

A cr = a -1 + lim / - — - — dl(q), 
tu-»i- y 1 — u>g 
o 

where, as earlier, measure / is the weak limit of sampling measures for N — > oo 

where A is an arbitrary Borel subset of [0, 1]. 

In |21] for closed networks similar results are obtained for more general dependence of intensities of the 
queue lengths. 

The trick with the graph extension allows get rid of other restriction: that for given cross the mean duration 
of red light is the same for all directions. It is necessary then, instead of vertex i, corresponding to this cross, 
introduce several vertices (i,d), where d enumerates possible directions on cross i. This poses some restrictions 
on the service times Tj^ in the new vertices, like 

n,d = n- 

d 

We limited ourselves to the problem, when in the system at least one jam appear. It is interesting to find 
number of jams and the mean number of cars standing in jams. 

Relation with real life This model is convenient because all parameters can be statistically estimated. 
Namely, statistical estimates of the parameters Pij,/ii look like (for example for constant /x^) 

j 

where A^,-(T) is the number of cars in the time interval [0,T], choosing direction j on the cross i. 

Practically interesting is the optimization of traffic lights, that can be achieved by choice of Ti^d and by 
traffic lights synchronization (that is totally ignored in any exponential model), and by changing matrix P by 
helping the choice of route. 
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